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We study the equations for the evolution of cosmological perturbations in / {TZ) and conclude that 
this modified gravity model can be expressed as a dark energy fluid at background and linearised 
perturbation order. By eliminating the extra scalar degree of freedom known to be present in such 
theories, we are able to characterise the evolution of the perturbations in the scalar sector in terms 
of equations of state for the entropy perturbation and anisotropic stress which are written in terms 
of the density and velocity perturbations of the dark energy fluid and those in the matter, or the 
metric perturbations. We also do the same in the much simpler vector and tensor sectors. In order to 
illustrate the simplicity of this formulation, we numerically evolve perturbations in a small number 
of cases. 


I. INTRODUCTION 

In the past few years there has been a growing realisation that dark energy [1] and modified gravity theories [2] 
models need to be confronted with observational data in a systematic way. This has generated interest in constructing 
frameworks and formalisms for comparing models, or classes of models, to data as opposed to testing individual models. 
There are number of approaches which have been developed to do this including the Effective Field Theory for dark 
energy [3-6], Parametrized Post Friedmann framework [7-10], and the Equation of State for perturbations (EoS) 
[11-13]. These very similar ideas correspond to parameterizations at the level of the perturbed action, perturbed 
gravitational field equations, and the perturbed dark energy fluid equations, respectively. In this paper we will 
concentrate on the EoS approach. 

Each of these parameterization schemes can be used in two different ways. The first is to construct arbitrary dark 
sector theories and the second is to map from a given model to the observationally combinations. In the latter of 
these —“model-mapping” — EoS approach prescribed how micro-physical degree of freedoms in the model combine 
to affect the evolution of quantities such as densities and velocity fields that are related to observables in terms of 
equations of state for the gauge invariant entropy and anisotropic stresses. Using the former approach preliminary 
constraints have been discussed in [14-16] based on presently available cosmological data including that from the 
Cosmic Microwave Background (CMB), weak leasing and redshift space distortions (RSDs). 

One of the most popular modified gravity theories is the f{R) class of models [17, 18]. The / iJZ) models of gravity 
are constructed by replacing the Ricci scalar in the Einstein-Hilbert action by an arbitrary function of the Ricci scalar 
/ (77.). Such models are well known to lead to an extra scalar degree of freedom and it has been shown, for example 
in [19-21], that observationally acceptable models can be constructed. 

In this paper we provide expressions for the equations of state for perturbations which completely characterize the 
linearized perturbations in / (77) modified gravity, including the scalar, vector, and tensor modes. Essentially we 
will model-map this theory into the EoS formalism. In doing this we show that the /(77) modification to General 
Relativity (GR) can be formulated as a dark energy fluid at background order —something which is well known — 
and also at first order in perturbations. As well as providing a physical interpretation and allowing these models to 
be included under the umbrella of the EoS formalism, we will see that writing the theory in this way can allow for a 
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very simple inclusion in codes, such as CAMS, used to calculate cosmological observables. 


II. BACKGROUND FIELD EQUATIONS 

The / {TV) models of gravity are characterized by the action 

S=\j {n + f{TV)} + S^, ( 2 . 1 ) 

where TZ is the Ricci scalar and S'„ is the action describing the standard matter fields. Natural units, c = h = Mpi = 1, 
are used throughout this paper. Varying the action (2.1) with respect to the space-time metric yields the field 
equations, 

( 2 . 2 ) 

where is the Einstein tensor and is the stress-energy tensor of the standard matter fields. All contributions 
due to / {TV) are packaged into the extra-term which we call the stress-energy tensor of the dark sector, explicitly 
formulated as 


\ - {Rfiu + g^zuGi - V^Vp) /k, (2.3) 

where is the Ricci tensor, and Direct calculation shows that is covariantly conserved, = 0, 

as is required by the conservation of the matter energy-momentum tensor, The background geometry is assumed 
to be isotropic and spatially flat, with a line element written as ds^ = —dt^ + a (t)^ Sijdx'^dxV where a {t) is the scale 
factor. Instead of the first and second order time derivative of the Hubble parameter, H, we use the dimensionless 
parameters^ 

h' 

= eH = -^, (2.4) 

where the prime denotes derivative with respect to d/d In a. The dark sector can be viewed as a fluid, with energy 
density Pde = and pressure The field equations (2.2) can be recast as 

f2„-|-Qde = l, -I- Wdef^de = fCn - 1, (2.5) 

where Q, = and w, = Pi/pi for i e {m, de}. From (2.3), the density and equation of state parameters of the f{TV) 
fluid are explicitly given by 

= -^ + (l-eH)/7^-4, (2-6a) 

Wde + 1 = — {‘^^nfn + (1 + ^H)fn — fn) ■ (2.6b) 

The case of a cosmological constant is recovered when f{Tl) = —2A, for which (2.2) reduces to the standard Einstein’s 
field equations. Note that (2.6a) is actually a second order differential equation for the function f{Tl) since f-ji = j' jTV. 
When the equation of state parameters, Wi’s, are taken to be constant, this equation can be integrated leading to the 
so-called designer /(T?.), see [22] for details. 


III. GAUGE INVARIANT FORMALISM FOR LINEAR PERTURBATIONS 


The dynamics of linear perturbations is written in Fourier space, in both the synchronous and conformal Newtonian 
gauges. Instead of the coordinate wavenumber that appears in the Fourier transform, k, a reduced dimensionless 
wavenumber will be used. 



(3.1) 


— ^ch). Furthermore, ch and eh are related through eh = — 2eh. 


^ With these notation the Ricci scalar reads 7^ = 
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so that K <C 1 and 3> 1 can be used identify the sub-(super)-horizon regimes. In the synchronous gauge, the non-zero 
metric perturbations are Sg^j = a^hij. In an orthonormal basis {k, l,m\ in fc-space, the spatial matrix /ly is further 
decomposed as hij = \h5ij + h\\crij + ■ Vij -I- ■ e^, where the notations contain the two vector (tensor) 

polarization states (the dot product is to be understood as a sum over the polarization states). Instead of h, we use the 
combination 6ry = /ly — h. The basis matrices are aij = kikj — ^6ij for the longitudinal traceless mode, = 2kf^jlj-^ 
and = 2k(^imj) for the vector modes and efj = efj = hlj — rhirhj for the tensor modes. In the conformal 

Newtonian gauge, the scalar modes are given by d^oo = —2a^'0 and 5gij = —2a^(j)Sij, while the tensor and vector 
modes remain the same in both gauges). An additional scalar degree of freedom arises at the perturbative level from 
the non-vanishing given by 


Ch 6i/2■ 


(3.2) 


This feature was pointed out in [23] and is a manifestation of the well-known connection between f(TZ) theories 
and non-minimally coupled scalar-tensor theories [24, 25]. Actually, constitutes the first non-trivial contribution 
of an arbitrary function of the Ricci scalar since a linear or affine f{Tl) can always be recast as standard GR with a 
rescaled Newton constant [22]. 

Our results are presented simultaneously in both the synchronous and conformal Newtonian gauges thanks to a 
new set of variables presented below: IT, A, Y, Z for the GR sector and y, y', x" for the f{TZ) sector. The introduction 
of this set of variables is motivated by the gauge transformation rules that are recalled in Appendix A. Quantities 
denoted with the subscript ‘S’ (‘C’) are evaluated in the synchronous (conformal Newtonian) gauge. 


Symbol 

Synchronous gauge 

Conformal Newtonian gauge 

T 

h'„ 

2K2 

0 

Y 

r -k e„T 

4’ 

Z 

V-T 

0 

X 

Z' + Y 

Z' + Y 

W 

A'-e„(A + F) 

X'-e^iX + Y) 

X 

X. + f^T 

Xc 

x' 

X's + (/^ - 

X'c - Ini’ 

x" 

Xs - CkXs + Ifn - Sch/k + (4eH - e^)fn)T 

Xc - EhXc - - 2(/k - eH/K)V’ 


Let us emphasize that x' and x" are not just the first and second derivatives of Xi but are degrees of freedom 
in their own right. Using these gauge invariant variables, the first order perturbation of the Ricci scalar reads 
STZ = + 4X — ^K^(Y — 2Z) — e^T). The fact that T appears explicitly in the expression of 5TZ indicates 

that this is not a gauge invariant quantity, and hence that x defined in (3.2) is also not gauge invariant. However, x 
defined in the table above is gauge invariant and it can be written in terms of the geometric perturbations as 

X=i{W + 4X-lK\Y-2Z)}, (3.4) 

which is valid in both gauges, when W,X,Y and Z are replaced by the corresponding expressions presented in the 
table. 

A generic stress-energy tensor, can be decomposed into 

SD^i, = {pS + SP) + {p + P) {u„6u^ + u^5uv) + SPS^u + P^^ui (3.5) 

where the density contrast is 5 = Spjp^ the Hubble flow is parametrized by Uy = (—1,0) in coordinate time, and 

buy = (0, dui) is the perturbed velocity field whose scalar mode is 0 = . 

Instead of b and 0, we make an extensive use of the dimensionless variables 

A = d-F 3 (1-f w) H6», 0 = 3 (1-k w) iJd. (3.6) 

In the same way as for the geometric perturbations, it is possible to form gauge invariant combinations of the perturbed 
fluid variables: 


Symbol Synchronous gauge 

0 e, + s{i + w)T 

bP bPs + P^T 


Conformal Newtonian gauge 

0c 

bPc 


(3.7) 
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The gauge invariant pressure perturbation, SP, is packaged into the gauge invariant entropy perturbation, 

wT= — -^(A-e). (3.8) 

p dp \ ) 

The anisotropic stress is the spatial traceless part of the stress-energy tensor. In the same way as the metric pertur¬ 
bation, it decomposes into one scalar, 11®, two vector, 11'', and two tensor modes, 11®’. Note that our B and H® differ 
from and a (anisotropic stress) defined in [26] by and {p P P)a = — |Pn®. 

The generic perturbed fluid equations which follow from the conservation of the stress-energy tensor, = 0, 

are 


A' — 3rcA — 2r(;n® -|- ^kCh© = 3 (1 -I- ic) A, 
0'-f3(^-w-f ieH)0-3^A-2wn®-3w;r = 3(l-fr(;)y, 

where 


- 1 

5k = 1 + T, -■ 

3e„ 


The field equations (2.2) expanded to linear order in perturbations, -|- yield 


- = n„A„ + flaeAa,, 

2A = n„0„ +11,3 0,3, 


fiT + 2A-|K2(r-z) = nj6Pjp^) + n,^{6P,jp^^), 

lK^{Y-Z) = + r!,3n;„nl, 

i/i"" + (i - len)h^' = n^w^ni + 

Ih^”+ + lK^h^ = n^w^ni + n,^w,^ni. 


(3.9a) 

(3.9b) 


(3.10) 


(3.11a) 

(3.11b) 

(3.11c) 

(3.11d) 

(3.11e) 

(3.11f) 


The first equation (3.11a) enables to write K^Z in terms of the A^’s, while the second equation (3.11b) consti¬ 
tutes the expression of the metric perturbation X in terms of the perturbed fluid variables 0i’s. The variables 
{W) A, y, Z, y, x', x", ft,'' ®'}, which are linear combinations of the metric perturbations and their time derivatives, 

are called the geometric perturbations. The variables ^A,, 0i, JP,, T,, 11®’''’®'| with i S {m, de}, which are linear 
combinations of the different projections of a perturbed stress-energy tensor, are called the perturbed fluid variables. 


IV. EQUATION OF STATE FOR PERTURBATIONS 


Equations of state for perturbations (EoS) constitute expressions for the entropy perturbation, T,,, and the 
anisotropic stresses, H®;''’®', that are constructed out of the perturbed metric degrees of freedom, the matter fluid 
variables, and the dark density and velocity divergence fields. Once these expressions are provided, the equations 
governing cosmological perturbations explicitly closes. Schematically, for the scalar sector in synchronous gauge, we 
are looking to obtain expressions of the form 


r,3 = r,3((5,3,5,3, ft', 77,..., (5„), n®^ = n®^((5,3, 5 , 3 , ft', 5 , • ■ •, <^m), 


(4.1) 


where the list of arguments shown is not exhaustive and can include derivatives, for example. Certain classes of 
equation of state have already been worked out (see [12, 13] for kinetic gravity braiding models, [27] for coupled 
Horndeski theories, [28] for generalised scalar-tensor theories and [29-31] for relativistic elastic and viscoelastic material 
models). 

The simplest way to understand this approach is in the vector and tensor sectors of the theory. If one assumes that 
there are no extra vector and tensor degrees of freedom (which is the case in f{TZ) theories), then the anisotropic 
stresses can only be functions of the metric variables. Focusing on tensor modes, the only tensor field available is the 
tensor mode of the metric perturbation, ft®, and its time derivatives (which we will limit to second order). The most 
general form of ll]}, would then be given by 

n® =rift^" + r 2 ft^' + r3ft^. 


(4.2) 
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where the {7i} are a set of dimensionless functions of space and time, that do not depend on the perturbed field 
quantities. Often one can deduce that these can be limited to just being functions of time only for specific theories. 
A similar expression could be written for the vector sector. 

In f{Tl) gravity, the expansion to first order in perturbations of the dark sector stress-energy tensor is 


+6 {VfjVvfn) - (n/T?,) - gf,„6 (Ofn) ■ 


(4.3) 


This allows us to isolate the perturbed fluid variables for the f{TZ) dark sector theory [32]. The tensor and vector 
projections of (4.3) readily constitute the EoS for and IIj^, 


V/ 


flaetaa.nl = -l{(3-e^)fn + fn}h 

flaetaa.nl = -ynh^'' - H (3 - + fn} 


(4.4a) 

(4.4b) 


As expected these are of the form (4.2) and the coefficients are just functions of time except for the explicit dependence 
on in the final term in the expression for H]]).. 

The scalar projections yield the following expressions: 


^^ae^ae — —9k^iiX + 

flaAe = x'-X-^MX, 

f^a»(^i^ae/Pa=) = gX" + (§ “ |eH)x'“ (l “ Ieh “ fK^) X 

-IfnW - 2{fn + |/4)^ + (E - Z) , 

flaetaa.n^ = -\K^x- Ifn^^Y - Z). 


(4.5a) 

(4.5b) 

(4.5c) 

(4.5d) 


From now on, the standard matter fluid will be assumed to have vanishing anisotropic stress and entropy perturbation, 
= r„ = 0 which is the case for a CDM fluid. When those are they are non-zero, the procedure presented below is 
easily generalized, with additional terms proportional to 11^ and r„. In the last equation (4.5d), 11^^ can be eliminated 
with (3.lid), providing the expression of Y in terms of Z and x, 


Y = Z-^X, 

valid for all K. Therefore, the dark sector anisotropic stress is simply 

^de'a'dellL = i+f^ ChX- 

Equation (4.5a), combined to (3.11a) and (3.11b), enables to write x in terms of the Aj’s and Oj’s, 


y =-_ 

Skch Skch 


^ {lIde(A,. - + fI„(A„ - ^0„)} . 


With (4.7, 4.8) one obtains the EoS for the anisotropic stress: 


w^deHl = 


3ffKCH 


K2 Ah„ - 


fiz 

2{l + fn) 


0d 


Oj_ 

tide l-t-Z-R 


^A„ - 


tim flZ 

fide 2(l + fTi) 


0. 




(4.6) 


(4.7) 


(4.8) 


(4.9) 


In order to deduce the entropy perturbation as an equation of state, one might begin with the expression for 
in (4.5c), and then eliminate x f’^d its time-derivatives, as in [12]. This involved differentiation of (4.5a) or (4.5b) 
to obtain y", and x'■ However, this strategy eventually leads to the perturbed fluid equation (3.9b) and hence a 
tautology, and therefore an alternative strategy is required. The starting point is the field equation (3.11c). On the 
right-hand-side of (3.11c), the pressure perturbations and 6P^ are replaced in favor of the Tj’s with (3.8). On the 
left-hand-side of (3.11c), W is replaced in terms of the geometric perturbations X,Y,Z and x with (3.4). As before, 
X is written in terms of the 0i’s with (3.11b). Furthermore, with (4.6) and (4.8), Y can be expressed with Z and 
the perturbed fluid variables A^s and 0i’s. After these replacements, the geometric perturbations only appear within 
K^Z which can be replaced with the A^s with (3.11a). Eventually, the EoS for the dark sector entropy perturbation 
is obtained as 


■f^derde 


A _ in 2{l+M-f^ 
Sgnin fj. 




c»- 


7^ 

Eh 2 / 7 ^ —/I 


SgRCH /4 


^^de CdeOde 


(4.10) 
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Figure 1: Evolution of some relevant f{7V) quantities that appear as coefficients in the EoS for the dark sector anisotropic 
stress (4.9) and entropy perturbation (4.10), for a designer f{TZ) that mimics A-CDM (wde = —1) with Bq = 1. When these 
are negative valued, their absolute value is plotted in order to accommodate the logarithmic scale. 


where 


c 


5k Eh ~ Eh 
35k 


dP, 

dPi 


Note that when the matter fluid is pressure-less and = — 1, then Eh = Ch and therefore 


A _ tgK-l 


c 


.<?K-1 

33k 


(4.11) 


(4.12) 


Equations (4.9) and (4.10), as well as (4.4a, 4.4b) for the vector and tensor sectors, are the main result of this 
paper. They constitute the EoS for perturbations in f{R) gravity expressed in a gauge invariant way. One could 
choose to express the EoS in terms of the dark sector perturbed fluid variables and the geometric perturbations X, Y 
and Z. Eor the entropy perturbation, this is achieved by replacing A„ and 0„ in (4.10) with (3.11a) and (3.11b). 
The dark sector EoS for perturbations are then expressed in a ‘self-consistent’ way, which does not depend explicitly 
on the perturbed fluid variables of the other fluid components 


= ^K2(y-z), 

= -^(A,,-0,j + |^i^(z-y)-3^(A + |K^z). 
The coefficients that play an important role in the EoS are either proportional to f^, and 

J-R. 

R =_£e_ 

«h(1-|-/7J.) ’ 


(4.13a) 

(4.13b) 


(4.14) 


or its inverse. Their evolution in the case of a designer f{TZ) with = —1 is plotted in figure 1. In order to illustrate 
some of the applications of the EoS formalism and the gauge invariant notations, we shall now describe the procedure 
for solving the linear perturbations in f(TZ) gravity. 


V. DYNAMICS OF LINEAR PERTURBATION IN /(7^) GRAVITY 

The dynamics of vector and tensor perturbations is straightforward to deduce and therefore we will only focus 
on the scalar sector. The dynamics of the scalar perturbations can be specified by writing the four perturbed fluid 
equations (3.9), plus one evolution equation for the geometric perturbation Z which follows from the definition of the 
gauge invariant notations (3.3), 


Z' = X -Y, 


(5.1) 
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Figure 2: Spectrum of the ratio ZjY (or —in the notation of [22]) for different values of the equation of state parameter 
when Bo = 1 (left) and different designer f{7V) scenarios parametrized by Bo and with Wde = —1 (right). On the x-axis, the 
wavenumber is written in units ‘/i/Mpc’, where h — 0.73 is the reduced Hubble constant. 


where X and Y are given in (3.11b, 4.6-4. 8 ) in terms of the perturbed fluid variables. In fact, the system of five 
differential equations is overdetermined: when K = 0, the field equation (3.11a) can be used to express in terms 
of A„,; and when K 7 ^ 0, the same equation gives Z in terms of the Aj’s. This approach is powerful and elegant: it 
provides an efficient way to solve the linear perturbation in f{R) gravity, and the phenomenology becomes transparent 
through the interpretation of the fluid variables. 

Before proceeding to this analysis, let us note that an essential feature of linear perturbation in f{TZ) gravity can 
be deduced from (4.6), when x is replaced by its expression (3.4) in terms of the geometric perturbations, 

BK^{Y - 2Z) = 3BW + 12BX + 3f^(r - Z), (5.2) 

where B is defined in (4.14). Since the geometric perturbations {W,X,Y,Z} shall remain bounded during their 
evolution, it appears that for K larger than B~^l'^ the ratio ZjY is driven to 1/2, as illustrated in figure 2 in the case 
of a designer /(7^). 

Let us consider the case of a matter fluid with = 11® = r„ = 0, along with a designer f{TZ) fluid that mimics 
A-CDM {Wde = —!)■ The function f{Tl) is determined by (2.6a). As shown in [22], the different solutions to (2.6a) 
can by parametrized by the single number B, defined in (4.14), evaluated today^ when Oq = 1 (analytical expressions 
for f(TZ) are available in some regimes [33]). For the numerical simulation we have chosen i?o = 1, and we have set the 
initial conditions for the perturbations at redshift z = 100, when B <C 1. At such high curvature, during the matter 
dominated era, the initial conditions for the perturbations must follow from the general relativistic expectation in 
order to be consistent with CMB observations. Hence, initially Aj, = ©a, = 0 and H„A„ = — |K^Z, = 2X, 

with X = Y = Z. As mentioned before, the five relevant dynamical equations are 

A'de = ~3Aa3 — gKCn^de ~ 2n®^, A'„ = —gKCnSm + 3A, 

©', = -sAde - Cneda - 2ni - SFde, ©L = -€„©„ +3y, (5.3) 

Z' = X-Y, 

where X and Y are replaced with (3.11b, 4.6-4. 8 ) in terms of the perturbed fluid variables, while H®^ and Fde are 
given in (4.9) and (4.10). Recall that in the conformal Newtonian gauge, Z = (p and Y = ip in our notations that 
follow [26], while Z — —4) and T = 'L in Song-Hu-Sawicki notation [22]. With this strategy we have successfully 
reproduced the results presented in figure 2 of [22], see figure 3. 

Alternatively one could have favored the geometric perturbations instead of the perturbed fluid variables. The 
way the equations have been written makes it straightforward to go from one picture to the other. Let us pick the 
conformal Newtonian gauge in order to illustrate this point. In the conformal Newtonian gauge, x' = x'c ~ (®®® 

section III). Therefore, equation (4.5b) combined to (4.6) and (5.1) yields a first order differential equation for the 


^ A subscript ‘0’ means that the quantity is evaluated today. 















e' 
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Figure 3: Evolution of the metric perturbations = —Z (left) and (right) for the ACDM expansion history 

(lOdo = —1, = 0.76). The dotted lines correspond to the ACDM scenario without f{TZ)- In this case the amplitudes of 

the perturbed potentials $ and "F- do not depend on the wavenumber. The black lines correspond to a f{TZ) scenario that 
mimics the cosmological constant {Wde = —1, Bo = 1). The amplitudes of the metric perturbations are now sensitive to the 
wavenumber Kq = ;^. Even the infrared limit, Ko = 0, represented by the dashed curves closest to the dotted lines, there is a 
disagreement with the ACDM predictions. 


geometric perturbation Y. The differential equation for X is simply provided by (3.4), where x is replaced with (4.6) 
and W with its definition, W = X' — CsiX + Y). Hence, the resulting set of equations to be solved is 

X' = (ch - 4)X - enH + i+fc+y Y iX, 

Y' = -X- 2j^Y - 0„' = -eH0„ + ST, (5-4) 

Z' = X-Y, 

valid in the conformal Newtonian gauge, where Z = (j),Y = tjj and X = (j)' + ijj. From here, the results in the 
synchronous gauge can be obtained using the gauge transformation rules. 

The last strategy would be to consider only the dark sector fluid variables and geometrical perturbations, eliminating 
0„ in the equation for Y in (5.4) with (3.11b) and taking the perturbed fluid equations for the dark sector. The 
resulting set of equations is 

+ (ch — 4)A + (2 — g^)enZ, 

z' = T^X + X-Z,^ 

x’ = ^de0de + + ‘^fnX + JtiZ, 

with and Td, given in (4.13). In much of the previous work, the linear perturbation in f{TZ) gravity have been 
solved with equations analogous to (5.5), with the difference that and 0de are replaced with A„ and 0„, thanks 
to the field equations (3.11), see [23]. 

We argue that the EoS approach (5.3) provides a clearer set of equations which can be solved and interpreted in an 
easier way. For instance, the stability of the metric perturbation in the high curvature regime, as discussed in [22], 
can be straightforwardly seen in (5.3). When H <C 1 (high curvature), from a quick look at the EoS (4.9, 4.10) we 
can see that 


0de — ~3Ad3 — eH0do ~ 2n®^ — 3rdj,, 

A'd, = -3Ad, - gKeH0de - 2ni, (5.5) 


n 


s 

de 


35k 


K^Ad 


Fd. = - 




(5.6) 


By considering the second derivative of the perturbed fluid variable Ade from (5.3) in the high curvature regime, 
one is left with a second order differential equation which can be recast as 

+ (3 - + I Ad, = |F(A„, 0„). (5.7) 

Therefore if i? < 0, the perturbed fluid variable Ad, does not converge toward the particular solution. Ad, = 
F(A„,,0„), but diverges exponentially, see [22] for a detailed discussion. 
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When a model of f{TZ) gravity is specified analytically, the EoS for and (4.9, 4.10) are easily obtained from 
a direct calculation of f-ji and f^. For instance, the Hu-Sawicki-Strarobinsky model [21, 34, 35], 


n 


, 2n 


fin) = -2A- 


commonly used as an alternative to the cosmological constant, leads to the following expressions for f-jz and 


(5.8) 


fn = 2 ri 


l-f:^ 


fm 




== 2 n( 2 n+l)- 


2n-l 

2n+l 




4n 


fin) 


l + 2 (^)''“ + 


2n 


n 


V 4n 


n n 


(5.9a) 

(5.9b) 


Another interesting case is when the dark energy fluid is dominating (n„ = 0, fide = !)■ Then the equation of state 
parameter reduces to 


1 + Wd, = fen- 


(5.10) 


This case is relevant not only for the present acceleration but also when studying inflationary scenarios based on fin) 
modifications to GR. In particular, when the slow-roll conditions are fulfilled, Ch <SC 1 and e(, ^ Ch, one gets 


dPde _ 
dpde 


and 


y _ 4 PK-l 
3 SK • 


(5.11) 


with = 2e| and Ch = (see footnote 1). With no other fluid than the fin) fluid, the perturbed field equa¬ 

tions (3.11) provide very simple relationship between the geometric perturbations and the perturbed fluid variables: 
— = Ade and 2X = For the most popular Starobinsky’s proposal [36] for primordial acceleration. 


fm 


7^2 

6 M 2 ’ 


(5.12) 


one finds f-ji = |, f^ = — | with Eh = , as well as = 1 + ^ 2^2 ■ So that the expressions for the anisotropic 

stress and the entropy perturbation become 

tCdeHl = (Ad, + eH0da) , (5.13a) 

^CdeTde = 3|^{(5K + i^)Ad. + (gK-l)0de}. (5.13b) 


These simple expressions can be plugged into the perturbed fluid equations (3.9) in order to solve the dynamics of 
linear perturbations during Starobinsky inflation, in the conformal Newtonian gauge or synchronous gauge. 


VI. DISCUSSION 

The EoS approach for dark sector perturbations has been discussed in details for (i) generalised k-essence theories, 
where the generic Lagrangian is Cdei4>,x)j with x = (ii) for theories in which the dark sector 

Lagrangian only contains the metric tensor C^^ig^^), see [29, 30]. In these two cases the gauge invariant equations of 
state were found to be 

= 0, Wd^Td, = (cg - Wd,)Ad^, (6.1) 

where Cg is the sound speed in the effective dark sector fluid, and 

U^d.n®^ = -|(Cs - ICde) |Ad. - 0de - 3(1 -F Wde)^! , WdeFd, = 0, (6.2) 


respectively. 

We have presented the EoS approach to cosmological perturbations in fin) gravity. After reviewing the formalism 
for describing the evolution of linear perturbation in fin) gravity, we have exhibited three equivalent ways to solve 
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their dynamics. In previous work, linear perturbation equations in f(TZ) gravity are solved in the geometric picture 
(5.4). Using the EoS approach (5.3-5.5) appears to have some advantages over the geometrical approach (5.4) because 
all the f(Jl) modification can be implemented in the dynamics by simply adding a new fluid species at the perturbed 
level, rather than modifying the whole set of equation for the geometrical variables. 

The main results of this paper are the equations of state for Tde and 11^, (4.4, 4.9, 4.10, 4.13). In these expressions 
the entropy perturbation and the anisotropic stresses are specified either in terms of the perturbed fluid variables of the 
dark sector and standard matter fluid, 11^^ = n®^(Ad 3 , 0^,, A„, 0„) and Tj,, = rd 3 (Ad,, 0^,, A„, 0„), or the perturbed 
fluid variables of the dark sector and the geometrical perturbations, 11®^ = nj^(y, Z) and Fa,, = rd,(Ad,, 0^,, X, Y, Z), 
thanks to the field equations (3.11). An important point is the extra degree of freedom, x, induced by a non-trivial 
f{Tl) modification to GR, is absent of these expressions. The elimination of this internal degree of freedom is the 
essence of the procedure. 

In order to illustrate the EoS formalism we have presented the EoS in the scalar sector for three different cases: 
(i) the designer f{TZ) in the high curvature regime (5.6), (ii) the analytical Hu-Sawicki-Starobinsky model for dark 
energy (5.9; to be plugged into 4.9, 4.10); and (iii) the Starobinsky proposal for inflation (5.13). 
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Appendix A: Gauge transformation rules 


Any expression written in the conformal Newtonian gauge can by translated into the synchronous gauge and vice 
versa. The relationship between both gauges can be seen as an infinitesimal coordinate transformation, with a specihc 
four vector characterizing the change of coordinate [26]. The time-like component of is 





(Al) 


It is related to T introduced in section III by T = HdP. The dot denotes derivative with respect to coordinate time. 
The transformation rules for the metric perturbations, from the conformal Newtonian gauge to the synchronous gauge 
are 


ijj = 

(j) = r] — Hd'^, 
^ + = rj — Hd'^, 

'(j) + + 2Hij} = ii-Hd°. 

For the fluid perturbations, the transformation rules are 

= (5, + 

P 

e, = + 

= 5Ps + P(f, 


n. 


= n. 


(A2) 


(A3) 


where the subscripts ‘c’ and ‘s’ hold for conformal Newtonian gauge and synchronous gauge respectively. For y and 
its time derivatives the transformation rules are obtained from its definition in terms of the first order perturbation 
of the Ricci scalar (3.2) 


Xc 

Xc - fnP 

Xc - fnP - ‘2fn'4’ 


Xs + fnd°, 
Xs + fnd°, 
Xs + fnd°- 


(A4) 


Writing these relations with T instead of d^ and the ‘prime’ derivative instead of the ‘dot’ derivative leads to the 
definition of the gauge invariant notations of section III. 
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